ANGULAR MOMENTUM TRANSPORT IN SOLAR-TYPE 
STARS: TESTING THE TIMESCALE FOR CORE-ENVELOPE 

COUPLING 

Pavel A. Denissenkov^'^, Marc Pinsonneault^, Donald M. Terndrup^, and Grant Newsham^ 

ABSTRACT 

We critically examine the constraints on internal angular momentum trans- 
port which can be inferred from the spin down of open cluster stars. The rotation 
distribution inferred from rotation velocities and periods are consistent for larger 
and more recent samples, but smaller samples of rotation periods appear biased 
relative to vsini studies. We therefore focus on whether the rotation period 
distributions observed in star forming regions can be evolved into the observed 
ones in the Pleiades, NGC2516, M34, M35, M37, and M50 with plausible as- 
sumptions about star-disk coupling and angular momentum loss from magnetized 
solar-like winds. Solid body models are consistent with the data for low mass fully 
convective stars but highly inconsistent for higher mass stars where the surface 
convection zone can decouple for angular momentum purposes from the radiative 
interior. The Tayler-Spruit magnetic angular momentum transport mechanism, 
commonly employed in models of high mass stars, predicts solid-body rotation 
on extremely short timescales and is therefore unlikely to operate in solar-type 
pre-MS and MS stars at the predicted rate. Models with core-envelope decou- 
pling can explain the spin down of 1.0 and 0.8 solar mass slow rotators with 
characteristic coupling timescales of 55 ± 25 Myr and 175 ± 25 Myr respectively. 
The upper envelope of the rotation distribution is more strongly coupled than 
the lower envelope of the rotation distribution, in accord with theoretical pre- 
dictions that the angular momentum transport timescale should be shorter for 
more rapidly rotating stars. Constraints imposed by the solar rotation curve are 
also discussed. We argue that neither hydrodynamic mechanisms nor our revised 
and less efficient prescription for the Tayler-Spruit dynamo can reproduce both 
spin down and the internal solar rotation profile by themselves. It is likely that 
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a successful model of angular momentum evolution will involve more than one 
mechanism. Further observational studies, especially of clusters younger than 
100 Myr, will provide important additional constraints on the internal rotation 
of stars and could firmly rule out or confirm the operation of major classes of 
theoretical mechanisms. 

Subject headings: stars: evolution — stars: interiors — Sun: rotation 



1. Introduction 



Rotation is an important attribute of the life of a star. When it is fast enough, rotation 
can trigger various (magneto-)hydrodynan iic instabili t ies that will drive mixing and angu- 
lar momentum transport in the star (e.g., IZahrull992l : ISpruitlll999l ). Unfortunately, except 
for the Sun, observations give information only about surface rotation of stars. Therefore, 
indirect methods have to be used to study rotation-driven transport processes in stellar 
interiors. When stars reach the main sequence (MS) their surface response to the torque 
from a magnetized wind depends on the timescale for internal angular momentum trans- 
port. If open clusters are treated as an evolutionary sequence, the assumptions of solid 
body (SB) and differential rotation (DR) lead to statistically distinguishable differences 
in the time evolution of the distributions o f cluster star rotation rates (see, for exa mple, 
Keppens. MacGregor. &: Charbonneaul Il995l : iKrishnamurthi et al.l 119971 : lAUainl Il998l ) . At 
late ages, models also have to be consistent with the strong coupling evident in the solar 
i nternal rotation profile, with nearly SB rotation in the r adiative core down to ~0.2i?Q 
( jTomczyk. Schou. fc Thompson! Il995l : ICouvidat et al.ll2003l ). 



The prior estimates of the timescale for core-envelope coupling in solar-type stars usually 
agreed on a short coupling time of order one Myr for the fastest rotators. However, they 
considerably disagreed on the coupling time fo r slowly rotating solar analogs for which the 



estimated va lues varied from 1 0-100 Myr (e.g., Keppens. MacGregor. &: Charbonneaulll995 



Allainlll998l ) to 0.5 Gyr (e.g., Ilrwin et al.l 120071 ). This discrepancy has been caused, on the 



one hand, by not carefully treating the selection effects and statistics and, on the other 
hand, by making particular assumptions about such things as the disk-locking time, initial 
rotation period distribution or the relative cluster ages which have since been updated. For 
example, the early work assu med that a Per was 50 Myr and the Pleiades was 70 Myr old 
( iMacGregor fc Brenneiill99ll ). In this work we apply rigorous methods of statistical analysis 
to the most recent large observational datasets on rotation periods of low-mass stars in open 
clusters and extensive computational modeling of their rotational evolution to constrain the 
timescale for core-envelope coupling in the slowest rotators. 



- 3 - 



In this paper we find that the timescale for core-envelope couphng depends on both 
rotation rate and mass, but it is significantly longer than would be predicted if SB rotation 
was enforced on a short timescale. A ngular momenturn transport by magnetic torques gen- 



erated by the Tayler-Spruit dynamo (ISpruitlll999l . |2002| ) does not pass this key test because 



it always enforces SB rotation in a solar-type star, no matter how slowly it rotates. This 
disagrees with the observational evidence that slow rotators in young clusters are most likely 
to possess DR (their cores rotate faster than their envelopes) rather than follow the P-age 
relations computed using SB rotation models. Furthe rmore, we show that even a revise d 



prescription for the Tayler-Spruit dynamo proposed by iDenissenkov &: Pinsonneaultl (120071 ). 
which reduces the effective magnetic viscosity by nearly two orders of magnitude, is still too 
efficient in redistributing angular momentum in the outer parts of the radiative core; this 
is inconsistent with the observed evolution of slowly rotating solar analogues in the P-age 
plane. The revised prescription comes to a better agreement with observations when it is 
supplemented by an additional angular momentum transport mechanism that should be able 
to operate in the inner core on a longer timescale. Such a mechanism is also needed to assist 
the revised prescription in shaping the solar SB rotation. 

In our work, three different theoretical models are employed to study the rotational 
evolution of solar-type stars: a simple double-zone model and two full stellar evolutionary 
models, one with a constant viscosity and the other with the effective magnetic viscosities 
provided by the original and revised prescriptions for the Tayler-Spruit dynamo. The models 
are described in detail and compared with each other in Section 2. Because the choice of 
basic model parameters is necessarily constrained by observations, we find it helpful to 
briefly introduce the samples of rotation period and v sin i data, that will be used for a more 
detailed analysis later in Section 3, already at the beginning of Section 2. In Section 3, we 
present the results of our computations of rotation period distributions for solar-type stars, 
as they evolve from the pre-main sequence (pre-MS) deuterium birth line to the solar age, 
and compare them with observations of rotation periods in open clusters. A brief discussion 
and our main conclusions are given in Section 4. 



2. Basic Theoretical Models and Observational Data 

For studying the angular momentum evolution of solar-type stars, it is important to 
specify correct initial and boundary conditions as well as to incorporate into the stellar 
evolution code three principal processes that are believed to govern changes of their surface 
rotation with time: disk-locking, internal transport of angular momentum, and angular 
momentum loss from the surface. In this section, we briefly describe common ingredients of 
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our rotating stellar evolutionary models, such as the input physics, the initial and boundary 
conditions, the law used for angular momentum loss, and our considered mechanisms for 
internal angular momentum transport. 

Because of the complex inter-relashionship between the physical processes involved in the 
rotational evolution of solar-type stars it is also important to use all available observational 
data that can constrain it. Fortunately, thanks to ongoing and planned ground and space 
based planet transit searches (such as the Deep MMT Transit Survey, the Monitor Project, 
KEPLER and COROT), extensive datasets of rotation periods for young and intermediate- 
age open cluster low-mass stars have been accumulating quickly during the last few years. 
The data that will be used in our paper to constrain the models are summarized in Table [TJ 
Their corresponding rotation periods are plotted in Fig. [T] (crosses) for illustration. 



2.1. Common Model Ingredients 



In our full rotating evolutionary computations we emp loy an upgraded version of the 
computer code used by iDenissenkov fc VandenBergl (120031 ). The most recent update is 
the adoption of Alan Irwin's improved EO^. In addition, the energy losses du e to neu- 
trino emission a re now calculated with the code distributed by lltoh et al.l (119961 ). We use 
OPAL opacities (IRogers fc Iglesiaslll992l ) for temperatures above ~ 10^ K, complemented by 
Alexander fc Ferguson! (119941 ) d ata for lower tempe ratures. Nuclear reaction rates are taken 
from the NACRE compilation (lAngulo et al.lll999l ). Gravitational settling is not included. 



We accept Zahn's concept of the rotation-induced anisotropic turbulence in stellar ra- 
diative zones, with horizontal com ponents of t he turbulent viscosity strongly dominating 
over those in the vertical direction (jZahrull992l ). which assumes that the horizontal turbu- 
lence has erased the latitudinal differential rotation. This allows to consider the angular 
velocity as a function of radius alone. We take into account small corrections to the stellar 
structure equations arisi ng from the distortion of equipot ential surfaces by such shellular 
rotati on (for details, see IDenissenkov &: VandenBergl |2003| ). We used the ICrevesse &: Noels 
( 119931 ) mixture of heavy elements. Our code has been calibrated to reproduce the solar 
luminosity Lq = 3.85 x 10^'^ ergs~^ and radius Rq = 6.96 x 10^° cm at the solar age of 
^0 = 4.57 Gyr. This procedure yields the initial hydrogen mass fraction X = 0.708 for the 
solar heavy-element mass fraction Z = 0.018, and a mixing length a of 1.75. 



^We use the EOS code that is made pubhcly available at 
GNU General Public License. 



http : / /f reeeos . sourcef orge . net/ 



under the 
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2.2. Initial and Boundary Conditions 



Our computations start on the pre-MS deuterium birth line of iPalla &: Stahlerl (Il99ll ). 
For stars with masses M = 1 Mq and 0.8 M©, the estimated birth line radii are 4.8 i?© and 
2.9 i?0, respectively. We assume that convective regions rotate as solid bodies at all times. 
Rotation at the birth line is specified by the initial rotation period Pq. In practice, we begin 
with a non-rotating fully convective model that is located above the birth line. This model 
is evolved down to the birth line where it is spun up to P = Pq, which yields the initial 
model for our rotating evolutionary computations. The initial period values are taken from 
the observed period distributions of solar-type stars in the youngest stellar clusters from our 
compiled data samples: the Orion Nebula cluster (hereafter, ONC or Orion), NGC2264, 
and NGC2362 (Tabled]). The Kolmogorov-Smirnov test gives high probabilities for these 
distributions to have been drawn from the same real distribution (Fig. [2]). Because stars still 
possess deep surface convection zones at these young ages, different ways of achieving the 
same surface rotation rate at the end of this stage will have similar final outcomes. 

During its approach to the zero-age main sequence (ZAMS), the protostar contracts. 
If its total angular momentum Jtot conserved on the pre-MS then the star would spin up 
as a result of the contraction. However, observations show that a large number of stars 
arriving at the ZAMS rotate slowly, as if their angular velocity rather than angular mo- 
mentum has remained constant (see Fig. 3 in iRebull. Wolff, fc Strom! |2004| ). The most 
plausible explanation is that the interaction of protostars with their accretion disks extracts 
angular momentum from the central object, reducing or preventing spin-up as t hey con- 
tract. This can occu r from magnetic coupling betwe en the protostar and disk (IKoenigl 



199ll : IShu et al.l 119941) or through an enhanced wind (jMatt &: Pudrita 120051) . As in most 



other works (e.g.. iKrishnamurthi et al.lll997l : iBouvier. Forestini. &: Allairull997l : lAllairull998 



Tinker. Pinsonneault. fc Terndrupl 20021 ). we model the interaction between the protostar 
and accretion disk in a simple way, assuming that during an initial time interval < t < Td 
the interaction maintains P(t) = Pq. The disk-locking time is considered as a free param- 
eter. 

In our stellar evolution code the transport of angular momentum is treated as a diffusion 
process described by 



d 



dMr 



dn 

'dWr 



where d/dt is a derivative taken at a constant M,., and is a viscosity whose physical nature 
has yet to be identified in real stars. Equation ([1]) needs two boundary conditions. A natural 
initial condition for it is that the fully convective birth line model rotates as a solid body 
with f2(0, Mr) = r2c(0) = 27r/Po, where VL^ = 27r/P is the angular velocity of star's convective 
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envelope. The inner SB-rotation boundary condition dfl/dMr = is applied at the surface 
of the innermost mass shell used in our computations, at ~ 0. To derive the outer 
boundary condition that has to be applied at the bottom of convective envelope = Mbce, 
we integrate equation ([T]) from Mr = to M^. = Mbce, taking into account that the time 
derivative of the angular momentum of radiative core is 



Jr. 



d 2 
dt 3 



A/bc 



Afbco ^ 
It 



r^n) dMr = Jtot - Je 



Here, Je = le^e is the angular momentum of convective envelope, Je being the envelope's 
moment of inertia. Finally, we obtain the following upper boundary condition: 



2 2 ■ 2 

Jtot ~ le^e ~ 2 '"bce^cMbco — — 



dn 



(2) 



in which Jtot is the rate of angular momentum loss from the stellar surface. 



2.3. Angular Momentum Loss and Internal Transport 



For the an gular momentum loss from the surface, we adopt the magnetized stellar wind 
prescription of iKrishnamurthi et al.l (119971 ): 



tot 



R/Rq 



M/M, 







min 



(3) 



Here, f2sat is the velocity at which the wi nd is saturated. The parameter Jlsat i s known to 



strongly depend on the stellar mass (e.g., lAndronov. Pinsonneault. fc Sillsll2003l ). In Fig. [31 



we have used upper 90th percentiles of the fl^ distributions for our cluster sample to adjust 
the appropriate values of f2sat for the three mass bins that we use in this study. For internal 
angular momentum transport w e assume that fast rotator s behave as SB rotators, a result 



both of theoretical calculations (IKrishnamurthi et al. 



20001 ) and indicated by prior spin down studies ( llrwin et al 



1997 



Sills. Pinsonneault. fc Terndrup 



20071 : there were prior claims by 



the Monitor group). The error bars to the percentiles were calculated using bootstrap simu- 
lations by generating 1000 sampling distributions for each cluster. Our adjusted r2sa.t values 
(shown in each panel in Fig. [3]) are close to those reported by lAndronov. Pinsonneault. fc Sills 
(120031 ). In particular, we have found that the same value of Qsat = 8Qq can be used for 
either of our "solar-type" mass bins, 0.7 < M/M© < 0.9 and 0.9 < M/M© < 1.1, to re- 
produce reasonably well their corresponding 90th percent iles. This value lies between the 
estimates Qsat = QAQq and f^sat = 10.5Qq obtained by lAndronov. Pinsonneault. fc Sills 
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(120031) for stars with M = 0.8 M© and M = 1.0 Mq, respectively. For the "fully convective" 
M < 0.4 Mq bin, we ha ve adiusted the parameter r2sa.t = 2 -5 fl(7) . It can be compared with the 
value of 1.8 Qq used by lAndronov. Pinsonneault. &: Sillj (120031 ) as a saturation threshold for 
their 0.4 Mq model star. Note that in all our computations we use stellar models with masses 
appropriate for the considered mass bins, namely 0.3 M©, 0.8 M0, and 1.0 M0. The param- 
eter in equation ([3]) is calibrated by requiring that our 0.8 and 1.0 models have 
n^ = nQ = 2.86 X 10"^ rad s'^ (P© = 25.4 days) at the solar age. For stars with M < 0.4 Mq, 
we use our 0.3 M© model and the solar calibrated value of ~ 3.19 x 10^''cm^gs~^ that 
has produced the solid curve in Fig. [3^. 



2.3.1. Constant Viscosity 

In this section, we consider angular momentum transport in rotating stellar evolution- 
ary models with a constant viscosity v{t, M^) = uq. Although the assumption of constant 
viscosity does not reveal the physical mechanism responsible for the transport of angular mo- 
mentum in solar-type stars, it nevertheless permits us to estimate both an instructive min- 
imum value (z/o)min that still results in the solar SB rotation and a maximum value (z/o)max 
such that the rotational evolution with any value in excess of it looks identical to that with 
^0 = (i^o)max- We have computed the evolution of a rotating 1 Mq star for two combinations 
of the disk-locking time (in Myr) and initial rotation period (in days): (rd,Po) = (6,8) and 
(2,3). For either of these combinations, the computations have been done for the same set 
of values of (],at = 8^0 and uq = 2.5 x 10^ 5 x 10^ 10^ 2.5 x 10^ 5 x 10^ 10^ 10^ and 
10^ cm^ s~^. Results are presented in Fig. H] with solid and dashed curves for the first and 
second combination of initial conditions, respectively. Models with uq < (z^o)min ~ 5 x 10^ 
cm^ s~^ have a residual DR at the solar age inconsistent with helioseismic data. On the 
other hand, the rotational evolution of models with uq > (z/o)max ~ 10^ cm^ s~^ is almost 
indistinguishable from one another. 



2.3.2. Effective Magnetic Viscosities from the Tayler- Spruit Dynamo 



SpruitI (I1999U2OO2I ) has elaborated upon the finding by iTaylerl (119731 ) that toroidal mag- 



netic fields frozen into plasma in a stellar radiative zone are always subject to a pinch-type 
instability. A release of magnetic energy by this instability causes concentric magnetic rings 
to slide sideways, mainly horizontally, along the equipotential surfaces and, to some extent, 
also along the radius. Magnetic induction makes it possible for the unstable radial displace- 
ment to produce a weak poloidal field at the expense of toroidal magnetic field B^. If the 
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radiative zone rotates differentially, the poloidal field can be stretched around the rotation 
axis into a new toroidal field that will again be subject to the Tayler instability. Spruit's 
original idea was that these consecutive poloidal/toroidal field generations might sustain each 
other under certain circumstances, thus forming a dynamo loop. The Tayler-Spruit dynamo 
could drive some mixing through unstable radial displacements with an effective diffusion co- 
efficient r]e, as well as some angular momentum transport by magnetic torques proportional 
to the product BrB^ with an effective viscosity u^. However, when considering circum- 
stances under which his proposed dynamo might work. Spruit did not take into account 
a reduction of the unstable long-scale horizontal displacement by the Coriolis force, there- 
fore overestimating rj^ and u^. The Coriolis force has properly been taken into account and 



the tr ansport coefficients have been revised correspondingly by iDenissenkov &: Pinsonneault 



(120071). Further model develo pme nts related to the Tayler-Spruit dyn amo were proposed by 
Zahn. Brun. fc MathisI J2007h and lRiidiger. Gellert. fc Schultzl J2009h . 



Whereas the original prescription for the Tayler-Spruit dynar no has been shown to pro - 
duce a solar rotation profile in agreement w ith helioseismic data (|Eggenberger et al.ll2005l ). 
the revised prescription fails to do so alone (IDenissenkov fc Pinsonneaultl 120071 ). One of the 
objectives of the present work is to see if these prescriptions can yield the evolution of surface 
rotation consistent with the rotation period data for solar-type stars in open clusters. 

Let us summarize the basic equations for the magnetic tran sport coefficients in the Spruit 



mechanism. For details, the reader is referred to the paper of IDenissenkov &: Pinsonneault 



( 120071 ). For Spruit's original prescription, it is convenient to present the effective magnetic 
diffusivity and viscosity in the following forms: 



rie = a Ky , and i'e = a K- 



'A* 2 

-y ■ 



(4) 



Here, q = \d\nQ/d\nr\ is the rotational shear, and y is a solution of the 4th order algebraic 
equation 



ay 



ay' + (3y~2 = 0, 



where 



a 



and (3 = 2 



(5) 



(6) 



are dimensionless coefficients. In equations (jlj) and (El), we have used standard notations for 
the thermal diffusivity 



K 



AacT^ 



(7) 
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where k and Cp represent the opacity and the specific heat at constant pressure, respectively, 
and for the T- and /i-component of the square of the Brunt- Vaisala (buoyancy) frequency 



A^l = ^(Vad- Vrad), and Nl = gip 



p 



(91n/i 



dr 



In the last expressions, g is the local gravity. Hp is the pressure scale height, Vad and 
Vrad are the adiabatic and radiative temperature gradients (logarithmic and with respect to 
pressure), and /i is the mean molecular weight. The quantities 6 = — {din p/d\nT)p ^ and 
ip = {d\np/d\nfi)prp are determined by the equation of state. In particular, for the perfect 
gas law 6 = (f = 1. 

When using the effective magnetic viscosity z/ = z/e in the angular momentum transport 
equation ([T]), it is important to remember that the Tayler-Spruit dynamo keeps operating 
only as long as the effective diffusivity r/e remains larger than the magnetic diffusivity r/mag- 
When rje approaches ?7mag the poloidal magnetic field decays through Ohmic dissipation faster 
than it is generated by the unstable radial displacements of magnetic rings. Although theory 
says nothing about the behavior of z/g in this limit, we can formally estimate a minimum 
value of the effective magnetic viscosity that is reached when r/e = r/mag- To do this, we 
replace rj^ with ?7jnag in the first of equations (jlj) and then solve equations (jl]-[6]) with respect 
to q and z/g- The effective magnetic viscosity approaches its minimum value of 

(z.e)n.in = 6.96 X loM — — - ( ^-—^ — , cm^ s-^ (8) 



when the shear is reduced to 

The normalizations in the last two equations, for which we have used values typical for the 
solar interior, assume that all quantities are expressed in cgs units. We have also introduced 
the reciprocal of the Roberts number e = rj^^g/K -C 1. 

The ratio P defined by the second of equations ([6]) remains small (/? -C 1) everywhere 
in the star until the age of ~ 30 Myr because nuclear reactions have not yet built up a 
sufficiently strong //-gradient in the stellar core. Given that by this age Qe approaches its 
maximum value that turns out to exceed ~3f20 ~ 10-^ rads-^ in most interesting cases 
(see Section [31 and panel a in Fig. [3]), it is obvious from equation ([H]) that (z/e)min will be 
of order 10*^ cm^ s-^ everywhere in our model star during the first tens Myr of its evolution. 
Furthermore, since (3 continues to remain very small outside the radius r ~O.2i?0, where 
^ to the solar age, we can expect that (z/e)min will keep values of order 10®- 10^ 
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in the outer part of the radiative core even at older ages. These expectations are confirmed 
by our detailed computations. Taking into account that equation ([8]) gives only a lower limit 
for the effective magnetic viscosity and the inner core with r < 0.2 Rq contributes less than 
10% to the total moment of inertia of the radiative core, we anticipate that Spruit's original 
prescription should always bring about the P-age relations characteristic of SB rotators, like 
those obtained with our constant viscosity model for uq > (t'o)max- However, we have to 
ensure that in our computations the viscosity declines abruptly as soon as the shear q is 
reduced below its critical value given by equation This requirement may leave some 
residual DR in the core of our final solar model. 

In our full rotating stellar evolutionary models we have used an approximate computa- 
tional method. It assumes that, as soon as a radiative core develops in a pre-MS star and its 
rotation profile begins to deviate from uniform one, the effective magnetic viscosity is large 
enough everywhere in the core to potentially restore its uniform rotation. Our test compu- 
tations confirm this. However, when the redistribution of angular momentum by magnetic 
torques has led to g ~ gmin, the Tayler-Spruit dynamo ceases to work, and its related vis- 
cosity should be replaced with the molecular one z/moi <^ (t'o)min- As we have noted, theory 
does not describe how this transition from {i'e)mm to u^oi occurs. When an increase of fie 
at t > rd caused by the residual pre-MS contraction of the star and its angular momentum 
conservation finally gives way to an Qe decrease due to the surface loss of angular momen- 
tum with the magnetized stellar wind, a much stronger DR tends to accumulate in the core. 
However, as soon as the shear exceeds its critical value qmin, a large viscosity of order (z/e)min 
will resume the redistribution of angular momentum by magnetic torques smoothing out the 
angular velocity gradient until q drops below gmin again. Following this quahtative picture, 
we put into equation ([1]) z/ = (z/e)min multiplied by an exponential factor that cancels z/ when 
q approaches gmin- 

Fig. O compares the results of our computations using this method for (ra, Pq) = 
(6, 8), (6, 30), and (2, 3) (solid curves) with results that we obtained for the same initial con- 
ditions but applying the constant viscosities z/q = 10^ and 10^ cm^ s~^ (dashed curves). We 
have fisat = SQq throughout. The comparison shows that, for Spruit's original prescription, 
the rotation period evolution of the Sun coincides with that of uq = 10^ cm^ s~^ > (z/o)max; 
i.e., it always corresponds to the SB rotation case. However, the solid curve in the bottom 
panel in Fig. [6] demonstrates that, unlike the constant viscosity model with uq > (z/o)max, our 
solar model computed using Spruit's original prescription does contain a small differentially 
rotating core, as we expected. Unfortunately , its p resence cannot be revealed by available 
hel ioseismic data (see , howe ver, iGarcia et al.l 120071 ). This result agrees with that reported 



by lEggenberger et al.l (120051 ) 
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For Spruit's revised prescription (jPenissenkov &: Pinsonneaultll2007l ). the effective mag- 
netic diffusivity and viscosity are given by the following equations: 

After the substitution of rj^ = r^mag into the first of these equations, we find that the revised 
Tayler-Spruit dynamo ceases to work when the shear approaches a critical value 



A comparison of coefficients in equations ([9]) and (jTTl) shows that in the second case the 
residual DR in the solar model is expected to be much stronger. A minimum value of z/g can 
be estimated from the second of equations ffTOl) in which we put rj^ = r^mag and q = gmin- As 
a result, we get 

where all numbers are again expressed in cgs units. 

It is instructive to compare the minimum values of the effective magnetic viscosity 



(z/e)§|j and (z^e)5n^ given by equations (IHl) and ( fT2l) . respectively. Their ratio is 

(-)i - 41,2 f^V" fi±iV'^ cj^r" (m"' i^X". (13) 



The last two equations show that, even though the revised value of (t'o)min is less than the 
original one, it is still large enough for us to employ the same computational method that we 
used to study the evolution of rotating solar-type stars with Spruit's original prescription. 

In Fig. [HI results of our apphcation of Spruit's original and revised prescriptions in full 
evolutionary computations of rotating solar models are compared with each other as well as 
with a rotating model with no internal transport of angular m omentum (solid, dashed, and 



dotted curve, respectively). We confirm the conclusion made by lDenissenkov fc Pinsonneault 
( 120071 ). who used a crude approximation 



Jr r 



to construct an f^-profile in their model of the present-day Sun, that the revised prescription 
for the Tayler-Spruit dynamo produces a solar model with a large rapidly rotating core. 
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in contradiction with helioseismic d ata (compare the da shed curve with the filled circles 
representing observational data from lCouvidat et al.ll2003l in bottom panel). Our new result 
obtained here is that, in spite of this, the revised prescription leads to the rotational evolution 
that has even shorter periods at ages older than ~ 30 Myr than the SB rotational evolution 
obtained with Spruit's original prescription (compare the dashed and solid curves in top 
panel). This behavior is explained as follows. As we anticipated, the values of (z/e)min given 
by equation (fT2|l indeed turned out to be large enough for the transport of angular momentum 
from the core to envelope to occur on a short "SB rotation" timescale like those we got in 
our constant viscosity models with z/q > (t'o)max ~ 10® cm^ s"^ However, because of an 
early build-up of a steep critical fi-profile (like that shown with dashed curve in bottom 
panel), below which the Tayler-Spruit dynamo ceases to work, the amount of the core's 
angular momentum available for transport to the envelope is diminished more and more as 
the evolution proceeds beyond an age of ~30 Myr. As a result, if we used the value of the 
stellar wind parameter = 2.8 x lO^^cm^gs calibrated for Spruit's original prescription 
then the P-age relation for the revised prescription would be nearly parallel to the solid 
curve in top panel but it would be located at longer periods, hence the solar rotation period 
would be overestimated. To match the solar surface rotation, we had to reduce to 
1.25 X 10'^''cm^gs for the revised prescription. This has shifted the P-age relation toward 
its location shown with the dashed curve. Note that we have also done test computations in 
which rje from the first of equations ffTOl) was used instead of r/mag to estimate z/g. We have 
not found noticeable differences with the results obtained using the approximate method. 



2.4. Double-Zone Model 



The two zone model 
employed by many others 


was originally proposed by 


M 


acGreeoi 


(1991 


) and it has since been 




VlacGreeor & Brennerl 


1991; 


Keppens. MacGreeor. & Charbonneau 


1995; 


Siess & Livio 


1997; 


AUair 


ill998; 


Irwin et al. 


200 


7). Angular momentum transport be- 



tween the radiative core and convective envelope is parameterized as follows. The core 
and envelope, with moments of inertia and rotation rates /c, flc, h and Qe respectively, 
are assumed to rotate as solid bodies. If Qc > then the maximum angular momen- 
tum that can be transfered from the core to the envelope (AJ)max is estimated as a differ- 
ence between the core's initial angular momentum Jc = Ic^c and the angular momentum 
<^c,eq = -^c^eq the corc will havc when its angular velocity becomes equal to that of the 
envelope. This final equilibrium velocity fioq is determined from the angular momentum 
conservation: Jc + Je = Jc,eq + >/e,eq, or Jc^c + -^o^e = (-^c + -^o)^eq- So, wc have 



(^'-^)max T I T ^e)- 
-'c ~r Je 
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As a free parameter, the double-zone model uses the core/envelope coupling time Tc that 
defines the rate (AJ)jnax/Tc with which the angular momentum is transfered from the core 
to the envelope. An advantage of this model is that one needs to know only how Jc, /e, Tbce? 
and Mbce are changing with time for a particular star (the last two quantities are required 
to take into account a displacement of the bottom of convective envelope during the angular 
momentum transfer). If this information is available (from full evolutionary computations) 
then the star's rotational evolution is obtained as a solution of the following system of ODEs: 

-TT = + o '■^bce^eMbce, (14) 

at Tc o 

^ = ^^^-|rL^^cMbce + jtot, (15) 
at Tc o 

where Jtot should be replaced by expression ([3]). 

Because of its simplicity the double-zone model can be used for Monte Carlo simulations 
and other statistical studies, like those conducted in next section. However, before doing this 
we want to relate the constant viscosity model and the double-zone model to one another 
through their parameters z/q and Tc. In Fig. [71 values of the coupling time are adjusted so 
that the double-zone model simulates some of the P-age relations that we computed with 
the constant viscosity model using the same initial conditions. We have established the 
following approximate correspondence between Uq (in cm^ s ^) and (in Myr): (z/q, r^) = 
(2.5 X 10^ 90), (5 X 10^ 40), (10^ 20), (2.5 x 10^ 7), and (10^ 1). This means that the 
double-zone models with the coupling time longer than ~ 40 Myr correspond to the constant 
viscosity models that, by the solar age, still possess residual DR inconsistent with helioseismic 
data. On the other hand, the double-zone models with Tc < 1 Myr reproduce the P-age 



relations for SB rotators (cf. lAUainl Il998l ): they are equivalent to the constant viscosity 
models with uq > (t'o)max ~ 10^ cm^ s~^ as well as to the models in which angular momentum 
is redistributed by magnetic torques generated by the Tayler-Spruit dynamo in its original 
prescription. 



3. Comparison with Observations 

From a theoretical standpoint, the rotational evolution of low-mass stars is a very com- 
plex process. It depends on a number of parameters, such as the initial rotation period Pq, 



^As the morphology of P-age relations is slightly different for these models, the adjusted values of Tc are 
uncertain within ^ 10%. 
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disk-locking time Td, angular velocity threshold for the magnetized wind saturation f2sat, and 
the rate of angular momentum redistribution expressed in terms of the core/envelope cou- 
pling time Tc, constant viscosity uq, or as a function of stellar structure and other parameters 
when a particular physical mechanism is chosen to describe it. Because of this complex- 
ity, the progress in this field is primarily driven by constantly accumulating and improving 
observational data and their statistical analyses. 

A number of physical mechanisms have been proposed to explain the solar SB rota- 
tion, such as the smoothing of DR by the back reaction of the Lorentz force emerging 
from the generation of a toroidal magnetic field by shear ing of a preexisting poloidal field 
(iMestel &: WeissI 119871 : ICharbonneau fc MacGregorl Il993l ) , angul ar momentum redistribu - 
tion by magnetic torques generated by the Tayler-Spruit dynamo (jEggenberger et al.ll2005l ). 
or by internal gravity wav es excited by turbulent eddies in the solar convective envelope 
(jCharbonnel fc Talorul2005l : see, however, iDenissenkov et al.ll2008l ). However, none of these 
mechanisms has been shown to agree or disagree with available rotation period data for 
solar-type stars in open clusters. In this paper, we subject the Tayler-Spruit dynamo to 
such an observational test. The other mechanisms will be tested in our forthcoming papers. 
We will start with Spruit's original prescription. As it always results in the P-age relations 
similar to those obtained with the double- zone model having a short coupling time of order 
1 Myr (compare Figs. [5] and [7]), we will use the latter as its substitute model. The main 
advantage of this replacement is that the double- zone model computations are very fast, 
therefore they can effectively be used to perform extensive parameter-space investigations. 

Our first step in reducing the number of free model parameters is the adjustment of 
the value of fisat (see Section [273|1 . For the SB rotating (tc = 1 Myr) double- zone model 
to reproduce as close as possible the upper 90th percentiles simultaneously for all of our 
compiled data for each of the three mass bins, we had to choose Qgat = 8 Qq for our 
two "solar mass" bins and Qsa,t = 2.5 fi© for the "fully convective" mass bin (upper solid 
curves in Fig. [3], these curves are also plotted in Fig. [1], with Qe being transformed back 
to P). This procedure assumes quite naturally that the most rapidly rotating stars in the 
samples evolve as SB rotators. In fact, it is impossible to construct a double- zone model 
with DR that would fit the 90th percentiles for the intermediate-age (~100 Myr old) clusters, 
those located immediately behind the maxima on the fig vs. age curves, without making 
unreasonable assumptions about its parameters. 

It turns out that the lower 10th percentiles for the stars with M < 0.4 Mq can also be 
approximated reasonably well with a SB rotation evolution curve (the lower and upper solid 
curves in Fig. [St and Fig. [It, respectively). This result is expected because a turbulent eddy 
viscosity in these fully convective stars should redistribute angular momentum very quickly 
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(e.g., see iTinker. Pinsonneault. fc Terndrupl |2002| ). Having said that, we should mention 
a recent evidence that casts some doubt on this simple interpret ation. The fi rst period 



i 



measurements for five stars with masses below 0.5 Mq in Praesepe (IScholz &: Eisloffell 120071 ) 
appear to disagree with the predicted model trends, but this sample is not of sufficient size 
to draw firm conclusions. 

On the contrary, it turns out to be very difficult for a SB rotation evolution curve to 
approach all the lower 10th percentiles for stars in the mass bins centered at 1.0 Mq and 
0.8 Mq (dashed curves in panels a and b in Fig. [3]), unless one takes a disk-locking time 
well in excess of 20 Myr or starts with a very slowly rotating star. If the double- zone model 
parameters were not constrained by observational data but could be chosen at our will then 
any mechanism of angular momentum transport that persistently produces and maintains 
SB rotation, in particular Spruit's original prescription, could easily be brought in agreement 
with the open cluster rotation period data for slowly rotating solar-type stars. Indeed, in this 
case one could simply choose appropriate combinations of the initial period and disk-locking 
time, one or both of which having to be sufficiently long, such that the double-zone model 
with SB rotation (tc = 1 Myr) and with those parameters applied would embrace all the 
periods for slow rotators, no matter how long they are (dashed and dot-dashed curves in 
panels a and b in Fig. [1]). 

However, observations do not allow such an arbitrary choice of parameters of a rotational 
evolution model for solar-type stars. In particular, the typical disk-locking time has been 
estimated to lie in a range between 2 and 10 Myr. This result is based on the measuring 



excess, that traces dust, or an H„ emission line width that traces accretion ( 


Hillenbrand 


2005 


Lvo & Lawson 


2005; 


Bouwman et al. 


2006; 


Javawardhana et al. 


2006; 


Damianov et al. 


2007 


1. A recent statistical analysis of available data on the pre-MS circumstellar disks has 



led iMamajekl (120091 ) to the conclusion that "the fraction of young stars with optically thick 
primordial disks and/or those which show spectroscopic evidence for accretion appears to 
approximately follow an exponential decay with characteristic time ~2.5 Myr". Taken at its 
face value, this means that, on average, only 10% of active disks around pre-MS solar-type 
stars are expected to survive by the a ge of 5.8 Myr. Howev e r, the dispersion of the disk life 
times can be quite large. For instance, ISicilia-Aguilar et al.l (120091 ) have found that ~50% of 
18 members of the t] Cham cluster studied by them still show signatures of a circumstellar 
disk by the age of 8 Myr. Therefore, in our models we will use the maximum disk-locking 
time (rd)inax = 20 Myr as a safe upper limit. 

With the maximum disk-locking time limited by the value of 10 Myr or even 20 Myr, 
the SB rotational evolution can explain the longest periods in the intermediate-age clus- 
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ters only if it starts with initial periods that are much longer than those measured in the 
youngest clusters (dashed and dot-dashed curves in panels a and b in Fig. [1]). However, the 
assumption that individual open clusters had their unique distributions of Pq in the past, 
with very different statistics, does not seem realistic because the same laws of physics had 
most likely shaped them during the star formation. This conclusion is supported by the fact 
that the period distributions for solar-type stars in three of our youngest clusters, ONC, 
NGC2264 and NGC2362, have high KS probabilities of having been drawn from the same 
real distribution (Fig. [2]). Therefore, we will assume that the period distributions of stars 
in clusters of different ages f{t, P) represent an evolutionary sequence that started with the 
same distribution /(0,Po)- Besides, we will consider the initial distribution of disk-locking 
times to be flat and random in the interval < < ('rd)max with the maximum value of 
(rd)max = 20 Myr. As a proxy for f{0,PQ), we will take period distributions for the three 
aforementioned young clusters. If the initial rotation of cluster stars was much slower than 
predicted by the ONC data then we should see clusters arriving on the MS with very low 
rotation rates, e.g. systems with ages below 100 Myr should already have a lot of slow rota- 
tors. Small samples from very young open clusters exhibit no such trend, but larger sample 
sizes are required. 



The assumptions that we have just made are not novel. For instance, llrwin et al.l (120071 ) 
evolved the observed rotation rates of low-mass stars in NGC 2362 forward in time trying 
to reproduce some statistics of Qe distributions in older clusters. Like us, they employed 
a double-zone mo del. To compare the results of their computations with observations. 



Irwin et al.l (120071 ) calculated the lower 25th and upper 90th percentiles of the observed 
distributions of for solar-type stars in NGC 2362, IC 2391/IC 2602, a Per, M34, and the 
Hyades. These statistics have been chosen to characterize the slowest and fastest rotators in 
the selected c lusters. Using our double-zone model with Tc = 1 Myr, we confirm the conclu- 



sion made by llrwin et al.l (120071 ) that the rotational evolution of the fastest rotators among 
solar-type stars can be simulated closely enough with short coupling times characteristic of 
SB rotators (upper dashed curves in Fig. [8]). On the contrary, much longer coupling times 
of order = 50-150 Myr and = 100-300 Myr for the mass bins 0.9 < M/Mq < 1.1 and 
0.7 < M/Mq < 0.9, respectively, as obtained in our computations (solid curves in the same 
figure), are required to explain the rotational evolution of the slowest rotator^. 

From a theoretical standpoint, the inability of the SB rotation model, on the one hand, 
and the ability of the DR model, on the other hand, to reproduce the location of the slowest 
rotators in the fie-age plot is caused by their, respectively, inappropriate and appropriate 



Irwin et al.l (|2007l ) reported much longer coupling times because they did not recalibrate the wind pa- 



rameter Kyj for the slowest rotators. 
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characteristics of mapping of /(O, Pq) to f{t,P) for long rotation periods. Apparently, the 
still unknown physical mechanism responsible for the internal angular momentum transport 
in solar- type stars must have a property, which is most likely related to fie, of making ro- 
tation of the radiative core and convective envelope more and more decoupled as one goes 
from a faster to slower rotating star. For example, the timescale for internal angular momen- 
tum transport by hydrodynamic mechanisms becomes shor ter as the rotation rate increases 



in this fashion (jPinsonneault. Kawaler. fc Demarqud Il990l ). Assuming that a quantitative 



characteristic describing this property changes continuously with Q^. (or P), this should re- 
sult in qualitative differences between the true Qe distributions for evolved clusters and those 
obtained with the double-zone model for different values of r^. These qualitative differences 
can be put on a scale by comparing the observed distributions of Qc with the modeled ones 
using the Kolmogorov-Smirnov test. This procedure is more rigorous than just matching the 
percentiles, because it uses statistical information encoded in the entire distribution rather 
than only in a part of it. 

Given the aforementioned robust results concerning the fastest solar-type rotators, the 
basic assumption of our following statistical analysis is that the rotational evolution of stars 
with sufficiently short initial rotation periods can be described by the double-zone model 
with Tc = 1 Myr, whereas that of slower rotating stars needs Tc ^ 1 Myr. In principle, 
we ought to introduce some monotonically increasing function Tc{Po) into our double- zone 
model that would produce a smooth transition from short-period SB rotators to stars with 
progressing DR that had longer initial periods. However, given the simplicity of the double- 
zone model, such approach looks over-complicated. Therefore, we have decided to employ a 
simple step function 

'^=™ = I r. » 1 Myr. if P„ > ft. 

where Pc is a critical period. 

The main objective of our statistical analysis of open-cluster rotation period data for 
solar-type stars is to get estimates of the combination of parameters (Pc, Tc) from f|T6|) that 
give the highest KS probabilities of the hypothesis that an observed Qe distribution for an 
open cluster of age t and our theoretical fl^ distribution computed for the same age t have 
been drawn from the same real distribution. It is assumed that the initial period distribution 
/(0,Po) is provided by the youngest clusters from our data compilation: ONC, NGC2362, 
or NGC2264 (Table [1]). To do the analysis, we have performed extensive double-zone model 
computations densely covering relevant regions (squares) of the ((rd)max, ^c) parameter space 
for a number of Pc values. The resulting Pks contours are plotted in Fig. (TD] (the mapping 
of ONC to M35 for the mass bin centered at M = 1.0 M©), Fig. [H] (NGC2362 to M50 
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for M = 1.0 Mq), Fig. [H (the mapping of NGC2362 to M34 for the mass bin centered at 
M = 0.8 Mq), and Fig. [T1(NGC 2362 to M50 for M = 0.8 M©). 

The Pks contour patterns revealed in Figs. [TUl [TT], and [T2] clearly show a decrease of 
the most probable value of the coupling time with an increase of the disk-locking time, as 
expected. Indeed, a star whose pre-MS rotation rate was being kept constant, as a result of 
its interaction with a circumstellar disk, for a longer time would arrive at the ZAMS having a 
smaller amount of angular momentum stored in the radiative core. Therefore, even a shorter 
coupling time will not lead to a fast rotation of its convective envelope simply because there 
is not much of angular momentum left in the core to be transported to the envelope. 

From Fig. [2] and Table [T] we can see that there is a large number of stars in the youngest 
clusters with periods Pq < 12 days. Yet the Pks contour patterns in Figs. [TDl [TT], and [T2] 
first become apparent, as Pc increases, and then get dissolved well before the critical period 
reaches the value of 12 days. This means that the double zone model cannot reproduce 
the fie distributions in evolved clusters under the assumption that all stars are SB rotators. 
There is always a group of stars with rotation of the radiative core decoupled from that of the 
convective envelope. For the mass bin centered at 1.0 Mq, a lower limit for the coupling time 
can be estimated as Tc > 30 Myr for (rd)max = 20 Myr, and Tc > 55 Myr for (rd)max = 10 
Myr (Fig. [T0|) . Fig. [11] gives nearly 20 Myr longer coupling times. Probable values of the 
parameter Pc that divides stars into SB rotators and objects possessing DR range from 2.5 
to 5.5 days for these two mappings. However, what is more important for us here is that, 
in all considered cases, we can completely rule out the pure SB rotation evolution (Pc = oo) 
as a solution. Consequently, Spruit's mechanism or any other physical mechanism that can 
only produce nearly SB rotation should be rejected as a prescription for rotational evolution 
of solar-type stars. The true model should have a parameter (presumably related to Qc) that 
allows it to tune up its rotational evolution so that the latter would resemble the rotational 
evolution of the double-zone model with the coupling time changing from = 1 Myr to 
Tc ~ 55 Myr. For stars with masses in the interval 0.7 < M/Mq < 0.9, we estimate much 
longer coupling times: Tc > 150 Myr for (rd)max = 20 Myr, and Tc > 200 Myr for (rd)max = 10 
Myr (Fig. [12] and Fig. [131) . This result can be anticipated from Fig. [8] where we compare the 
lower 10th percentiles with our DR double-zone model computations. 

As an additional statistical exercise aiming to demonstrate the failure of SB rotation 
to serve as a unique solution for all solar-type stars, we have evolved the observed period 
distributions for M35 and M50 backward in time to zero age, assuming the short coupling 
time of 1 Myr for all stars, to see how their initial period distributions might resemble those 
of Orion and NGC 2264, respectively. To solve this inverse problem, we have first projected 
(using our double-zone model with Tc = 1 Myr) a sufficiently wide rectangular domain of 
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the (logPo; logTd) parameter space into the logP space at the ages of M35 and M50. The 
projections have turned out to look hke perfect planes. Contour lines for the M 35 projection 
are plotted in upper panel in Fig. O It is seen that the backward solution is not unique. In 
fact, every value of logP(t) at an older age t is mapped into a diagonal line segment 

logPo ~ alogTd + 6 (17) 

back at t = 0. To assign a unique value of the initial period for log P(t), we have decided to 
randomly select from an interval restricted by the end points of the diagonal (fTTIl . Thus 
computed "initial" rotation periods for the M 35 and M 50 stars are compared with the Orion 
and NGC2264 period distributions in Fig. M (middle and bottom panels). 

We see that the hypothesis that all solar-type stars, independently of their rotation 
periods, evolved like SB rotators would require that M35 and M50 had initially contained 
much larger fractions of slow rotators compared to Orion and NGC2264. Because this 
is unlikely to be the case, we should reject any prescription for the internal transport of 
angular momentum that enforces SB rotation in a solar-type star without respect to how 
slowly it rotates. In particular, the original Tayler-Spruit dynamo cannot be considered as 
a relevant mechanism for angular momentum redistribution in solar-type stars because it 
always enforces nearly SB rotation in both fastest and slowest rotators. 



4. Discussion and Conclusions 

One of the main conclusions made in this paper is that the observed rotation period 
distributions of low-mass stars in open clusters do not seem to support the hypothesis that 
all solar-type stars evolve as SB rotators. Instead, statistical analyses of these data show that 
only the fastest rotators among solar-type stars can be considered to possess SB rotation 
during their entire evolution, whereas their slowly rotating counterparts are most likely to 
manifest internal DR between t he ages of "^ 3 My r and several hundred Myr. This is not a 



new result though because, e.g., Ilrwin et al.l (120071 ) came to a similar conclusion. A novelty 



of our work is that we have used the entire period distributions of solar-type stars in a 
number of open clusters of different ages, not just some of their statistics, and extensive 
Monte-Carlo simulations to put this conclusion on a more rigorous quantitative basis. This 
seems to be quite a reasonable approach to the solution of the problem, given that it has 
a large number of poorly constrained parameters. In particular, we have found that a star 
with M = 1.0 ± 0.1 M0 that starts its rotational evolution with a period Pq > 2-4 days 
should have rotation of its convective envelope and radiative core coupled on a timescale 
of order Tc = 55 ± 25 Myr, where the systematic uncertainty of this estimate takes into 
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account the anticipated decrease of the couphng time with an increase of the disk-locking 
time. For a shghtly less massive star with M = 0.8 ± 0.1 Mq, the coupling time increases 
to Tc = 175 ± 25 Myr. Given that the initial period distributions, those for the youngest 
open clusters, have rather densely occupied bins up to Pq ~ 12 days (Fig. [2]), it turns out 
that quite large fractions of solar-type stars (up to 50%) should go through a phase of DR 
evolution. 

It is important to note that solar-type stars in open clusters older than a few hundred 
Myr are not suitable for a type of statistical analysis employed by us. The problem with 
the older clusters is that rotation periods of solar-type stars in them have already converged 
too close to each other, all aiming eventually to approach the solar rotation. Therefore, 
period distributions for the older clusters do not allow to make an unambiguous conclusion 
about the rotational evolution of solar type stars. To illustrate this, we have evolved the 
NGC2264 period distribution of stars with M = 1.0 ± 0.1 M0 to a period distribution at 
t = 550 Myr corresponding to the age of M37 (Fig. [T4|) . We see that at this old age the 
period distribution mapping admits two branches of solutions, one with DR and another 
with SB rotation. This bimodality is caused by the very narrow range of the mapping at old 
ages which finally degenerates into a point at the solar age. 

The second novelty of our investigation is that we have related the coupling time from 
the double-zone model to its corresponding value of the constant viscosity from our full 
stellar evolutionary model (Fig. [7]). In particular, the minimum coupling time of 30 Myr 
obtained in our analysis of DR of solar-type stars roughly corresponds to uq = 7.5 x 10^ 
cm^s~^, whereas the longest coupling time of ~80 Myr for stars with M = 1.0 ± 0.1 
implies that in some of them the internal transport of angular momentum can be as slow as 
that modeled with a viscosity z/q ~ 3 x 10^ cm^ s^^. 

Finally, we have shown that original Spruit's prescription always enforces SB rotation 
in a solar-type star, no matter how slowly it rotates. Therefore, it cannot be accepted as a 
physical mechanism for the internal angular momentum redistribution in radiative cores of 
solar-type stars many of which are indirectly proved to possess a degree of DR. 



We have also found that the revised prescription for the Tayler-Spruit dynamo (IDenissenkov fc Pinsonn 



20071 ) results in the rotational evolution superficially resembling that obtained using the orig- 
inal prescription, although the former leaves a large differentially rotating radiative core in 
the present-day solar model (Fig. [6]). Hence, the revised prescription appears to be in conflict 
both with the rotation period data for solar-type stars in open clusters (like the original pre- 
scription) and with the helioseismic data (unlike the original prescription). It is worth seeing 
if other possible angular momentum transport mechanisms can assist the revised prescrip- 
tion in bringing the rotating solar model closer to observations. To test this idea, we have 
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supplemented the effective magnetic viscosity f|T2|) with viscosities arising from the secular 
shear instability (denoted by the subscript "ss" below), Goldreich-Schubert-Fricke instability 
(GSF), as well as with the molecular viscosity (mol) and a viscosity = \rU{r)\/b that 
approximately describes the transport of angular momentum by meridional circulation as a 
diffusion process: 

^ = i'^c)hn + ^ss + Z^GSF + l^mc + ^'mol- (18) 

To calculate the quantities Ugs, ^gsf, and the meridional circul ation velocity U(r), we have 



used t he corresponding equations from Appendix to the paper of lChaname. Pinsonneault. fc Terndrup 



(120051 ) ■ Note that their expression for U{r) neglects terms depending on the /x-gradient and 
its derivatives. Therefore, the meridional circulation in this approximation is expected to 
penetrate deeper into the r adiative core than in th e case of its fully consistent implementa- 



tion originally proposed by lMaeder fc ZahnI (119981 ) that had later bee n applied to st u dy th e 



rotational mixing in solar-metallicity MS stars with M > 1.35 Mq by lPalacios et al.l (120031 ). 



In Fig. [ini the dotted curves represent results of our computations with the combined 
viscosity (ITS]) substituted into equation ([T]) . Note that a slightly increased amount of angular 
momentum that can be transported from the core to the envelope on a longer hydrodynamical 
timescale has changed our results in the right directions: first, after the wind constant K„ 
is properly recalibrated, the P-age relation has shifted toward longer periods, and, second, 
the degree of DR in the core has slightly been reduced. To find out if these changes will 
further grow in the same directions when the efficiency of supplementary angular momentum 
transport mechanisms increases, we have multiplied z/gsf, the dominating viscosity in the 
core, by a factor of 10. Results of this artificial viscosity enhancement are plotted in Fig. [15] 
with dashed curves. We see that both the P-age relation and the core rotation profile have 
indeed continued to change in the right directions. This exercise shows that the revised 
prescription cannot be rejected as easily as the original one. There remains a possibility 
that, when being assisted by other transport processes, it may reproduce the observational 
data yet. 

Our main conclusion is critically based on a comparison of theoretical predictions and 
observations of rotation periods for the slowest rotators in the intermediate-age open clus- 
ters. Therefore, if the observational data were biased toward the longest periods that would 
undermine confidence in our results. The real situation turns out to be opposite. Observa- 
tions for some of our used open clusters (for those with smaller data samples) are actually 
biased toward the shortest periods, as is expected from the photometric period measurement 
procedure that needs longer observational times to accumulate data sufficient for extracting 
longer periods. The following example illustrates this bia s. Filled circles in upper panels 



of Fig. [16] represent the Pleiades v sin i unbiased data from lAndronov. Pinsonneault. fc Sills 
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( 120031 ) ■ For comparison, open circles in the upper left panel show the Pleiades data used in 
our work that do appear to be biased toward lower rotation periods. On the other hand, 
open circles in the upper right panel are the M 35 data from iMeibom. Mathieu &: Stassun 



( 120091 ) ■ Their original periods have been transformed into vsmi values using a randomly 
generated angle < i < 7r/2 and R = Rq. The lower right panel shows that the unbiased 
data for the two clusters of similar age look alike (Pks = 0.112, the medians are 6.8 and 8.5, 
the first and third quartiles are 5.4 and 13 for the Pleiades, and 4.7 and 18 for M35). On 
the contrary, the biased Pleiades data in the lower left panel have very different median and 
third quartile values, 15 and 39, respectively. 

To sum up, our main conclusions can concisely be formulated as follows. Whereas the 
period distributions for the fastest rotators among solar-type stars in open clusters are very 
well reproduced assuming their SB rotation, the period distributions for the slowest rotators 
are better described by stellar models with P R in their radiativ e cores. This conclusion is in 



agreement with previous results reported by llrwin et al.l (120071 ). Our new result is that the 



original prescription for the Tayler-Spruit dynamo always enforces SB rotation in a solar- 
type star even if the star is a slow rotator. Therefore, this angular momentum transport 
mechanism is unlikely to operate in solar-type early MS stars. The revised prescription for 
the Tayler-Spruit dynamo cannot explain the observed period distributions either. Besides, it 
leaves a large rapidly rotating radiative core in the present-day solar model. To be consistent 
with observations, the revised prescription needs to be assisted by other angular momentum 
transport mechanisms that must be able to penetrate into the deep core and operate on a 
longer timescale. 



We are grateful to Luisa RebuU for providing us with the latest data on rotation periods 
of open cluster low-mass stars. We acknowledge support from the NASA grant NNG05 
GG20G. 
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Fig. 1. — Rotation period data (crosses) for the three mass bins used in this paper (see 
Tabled]). Curves represent the SB rotational evolution computed with different values of 
the initial period and disk-locking time using the double-zone model and common model 
ingredients that are described and discussed in text later. 
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Fig. 2. — Panels a and c: Period distributions for solar-type stars in the Orion cluster 
(shaded histograms) are compared to those in NGC 2264 and NGC 2362 (thick sohd curves). 
The Kolmogorov-Smirnov test gives high probabilities (Pks ^ 0.05) that the corresponding 
pairs of distributions have been drawn from the same real distribution. Panels b and d: 
Cumulative distribution functions for the two pairs of clusters. 
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Fig. 3. — Upper 90th and lower 10th percentiles of f2e = li^jP distributions (crosses) for 
all of our compiled data samples (Table [1]). The 90th percentiles are used to adjust the 
parameter f2sat (shown in each panel) by fitting them to the SB rotation evolution (solid 
curves). The 10th percentiles for the fully convective stars are also fitted well by a SB 
rotation evolution curve (panel c). The percentile vertical errorbars were evaluated through 
bootstrap simulations. 
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Fig. 4. — Rotation period evolution of the Sun computed using the constant viscosity 
model for the initial conditions (rd, Pq) — (6 Myr, 8 days) (sohd curves), and (ra, Pq) — 
(2 Myr, 3 days) (dashed curves), and for the values of i/q = 2.5 x 10^, 5 x 10^, 10^, 2.5 x 
10^, 5 X 10^, 10^, 10^, and 10^ cm^s~^ (from upper to lower curve for either combination 
of initial conditions) . 
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Fig. 5. — Rotation period evolution computed for the initial conditions (rd, Pq) = 
(6 Myr, 8 days), (6 Myr, 30 days), and (2 Myr, 3 days) using the original prescription for the 

Tayler-Spruit dynamo (solid curves) and the constant viscosity model with z/q = 10^, and 
10^ cm^ s~^ (dashed curves). The dashed curves for vq = 10^ cm^ s~^ almost coincide with 
the solid curves. All computations have been performed with M = 1 Mq and Ogat = SO©. 
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Fig. 6. — Rotation period evolution (upper panel) and final rotation profiles in the solar 
models (bottom panel) computed using Spruit's original (solid curves) and revised prescrip- 
tion (dashed curves). Dotted curves represent a model with no internal ang ular momentum 



transport. Filled circles in the bottom panel are helioseismic data from ICouvidat et al. 
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Fig. 7. — Comparison of rotation period evolution of the Sun computed using the con- 
stant viscosity (sohd curves) and double-zone model (dashed curves). The corresponding 
pairs of the viscosity and core/envelope coupling time adjusted for the models' P-age re- 
lations to resemble one another as closely as possible are {uq, Tc) = (2.5 x 10^, 90), (5 x 
10^ 40), (10^ 20), (2.5 X 10^ 7), and (10^ 1), where i/q is given in cm^s-\ and Tc in Myr. 
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Fig. 8. — Upper 90th and lower 10th percentiles (crosses) of Qc distributions for all of our 
compiled data samples of open cluster solar-type stars are compared with l^c-age relations 
computed using the double-zone model. Dashed curves correspond to Tc = 1 Myr (SB rota- 
tion), solid curves (from upper to lower) to DR with — 50, 100 and 150 Myr (panel a), 
and Tc = 100, 200 and 300 Myr (panel b). 
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Fig. 9. — Upper panel — contour lines in a domain of initial conditions [logPo, logTd) that 
are projected to specified values of logP (shown on the lines) at the age of M35 (150 Myr) 
for SB rotational evolution computed using the double-zone model with Tc = 1 Myr. These 
lines are used to restore the initial period distribution for solar-type stars in M35 (shaded 
histogram in middle panel) and to compare it with the Orion period distribution (thick solid 
curve). Bottom panel shows the corresponding comparison between M50 and NGC2264. 
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Fig. 10. — Contours of the Kolmogorov-Smirnov probability of the hypothesis that the 
observed and modeled period distributions of M35 stars with M = 1.0 ± 0.1 Mq have been 
drawn from the same real distribution. The initial periods are taken from the ONC cluster. 
The double-zone model with the dependence f lT^ of on Pq has been employed. Simulations 
have been done for the shown intervals of the parameters Tc and (rd)jnax5 and for the values 
of Pc specified in parenthesis atop each panel. 
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Fig. 11. — Contours of the Kolmogorov-Smirnov probability of the hypothesis that the 
observed and modeled period distributions of M50 stars with M = 1.0 ± 0.1 Mq have been 
drawn from the same real distribution. The initial periods are taken from the NGC2362 
cluster. The double-zone model with the dependence ( TTB]) of r^. on Pq has been employed. 
Simulations have been done for the shown intervals of the parameters and (Td)niax; and 
for the values of Pc specified in parenthesis atop each panel. 
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Fig. 12. — Contours of the Kolmogorov-Smirnov probability of the hypothesis that the 
observed and modeled period distributions of M34 stars with M = 0.8 ±0.1 Mq have been 
drawn from the same real distribution. The initial periods are taken from the NGC2362 
cluster. The double-zone model with the dependence ( IT6l) of Tc on Pq has been employed. 
Simulations have been done for the shown intervals of the parameters and (rd)niax; and 
for the values of Pc specified in parenthesis atop each panel. 
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Fig. 13. — Contours of the Kolmogorov-Smirnov probability of the hypothesis that the 
observed and modeled period distributions of M50 stars with M = 0.8 ±0.1 Mq have been 
drawn from the same real distribution. The initial periods are taken from the NGC2362 
cluster. The double-zone model with the dependence ( IT6l) of Tc on Pq has been employed. 
Simulations have been done for the shown intervals of the parameters and (rd)niax; and 
for the values of Pc specified in parenthesis atop each panel. 
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Fig. 14. — Contours of the Kolmogorov-Smirnov probability of the hypothesis that the 
observed and modeled period distributions of M37 stars with M = 1.0 ±0.1 Mq have been 
drawn from the same real distribution. The initial periods are taken from the NGC2264 
cluster. The double-zone model with the dependence ( IT6l) of Tc on Pq has been employed. 
Simulations have been done for the shown intervals of the parameters Tc and (rd)max, and 
for the values of Pc specified in parenthesis atop each panel. Because of the convergence to 
the solar rotation, there is an ambiguity in the rotational evolution mode (SB or DR) for 
this 550 Myr old cluster. 
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Fig. 15. — Rotation period evolution and final rotation profiles in the solar models computed 
using Spruit's revised prescription (solid curves) and the combined viscosity (ITSil (dotted 
curves). Dashed curves represent results obtained with the combined viscosity in which the 
term z/qsf has been multiplied by a factor of 10. For comparison, dot-dashed curve in upper 
panel shows the rotational evolution computed using the double-zone model with = 100 
Myr (a model with DR). 
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Fig. 16. — Filled circles in upper 



panels are the Pleiades v sin i unbiased data from 



Andronov. Pinsonneault. fc Sillj (120031 ) . Open circles in the upper left panel are the Pleiades 



data used in this paper that appear to be biase d toward lower rotation periods. Op en circles 
in the upper right panel are the M 35 data from lMeibom. Mathieu fc StassunI (120091 ). Periods 
have been transformed into v sini values using a randomly generated angle < i < 7r/2 and 
R = Rq. The lower right panel shows that the unbiased data for the two clusters of similar 
age look alike (Pks = 0.112, the medians are 6.8 and 8.5, the first and third quartiles are 
5.4 and 13 for the Pleiades, and 4.7 and 18 for M35). On the contrary, the biased Pleiades 
data in the lower left panel have very different median and third quartile values, 15 and 39, 
respectively. Vertical dotted lines in upper panels show the range of 0.7 < M/Mq < 1.1. 



Table 1. Rotation Period Data 








Number of Stars 


in the Sample 




Cluster 


Age (Myr) 


M/Mq < 0.4 0.7 < M/Mq 


< 0.9 0.9 < M/Mq < 1.1 


Source 


ONC 


2 


32 


3 


53 


RebuU. Wolff. & Strom f2004) 


NGC2264 


4 


13 


12 


63 


Rebuff. Wolff. & Strom r2004) 


Lupus 


4 





2 


12 


Rebuff. Wolff. & Strom r2004) 


NGC 2362 


5 


70 


46 


44 


Irwin et al. ('2008b) 


Tau-Aur 


6 





2 


23 


Rebuff. Wolff. & Strom r2004) 


NGC 2547 


30 


94 


26 





Irwin et al. ('2008a) 


IC2391 & IC2602 


40 





8 


28 


Rebuff. Wolff. & Strom (2004) 


a Per 


70 





8 


28 


Rebuff. Wolff. & Strom r2004) 


Pleiades 


115 





20 


19 


Rebuff. Wolff. & Strom (2004) 


Pleiades (wsini) 


115 


39 


50 


51 


Andronov. Pinsonneault. & Sills f2003) 


M50 


130 


163 


252 


62 


Irwin et al. f2009) 


NGC 2516 


180 


183 








Irwin et al. C2008b) 


M35 


150 





154 


111 


Meibom. Mathieu & Stassun f2009) 


M34 


220 


10 


35 


8 


Irwin et al. f2006) 


M37 


550 





197 


128 


Hartman et al. (2009) 


Hyades 


600 





8 


17 


Radick et al. fl987): Prosser et al. (1995) 



